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ABSTRACT 

The  scanning  tunneling  microscope  (STM)  and  the  atomic  force  microscope 
(AFM)  were  both  used  to  study  the  evolution  of  the  morphology  of  graphite  (0001) 
surfaces  etched  with  5  ke\'  Ai'^  ions  as  a  function  of  both  ion  beam  flux  and 
fluence.  STM  and  AFM  topographs  collected  after  ion  beam  exposure  show- 
surface  features  that  increase  in  average  size  as  a  function  of  ion  beam  fluence, 
but  which  have  a  non-monotonic  dependence  on  ion  beam  flux.  Furthermore,  the 
surface  roughness  does  not  follow  a  square  root  dependence  with  respect  to  the  ion 
fluence,  as  expected  for  a  stochastic  process.  One  explanation  for  this  behavior  is 
that  the  sputter  yield  for  graphite  is  a  function  of  the  ion  beam  fluence. 


Ion  beam  etching  (IBE)  is  a  widely  used  technique  for  surface  preparation 
and  analysis.!  Nonequillibrium  surfaces  produced  by  this  technique  have 
topographies  that  are  known  to  depend  on  the  sputtering  conditions. ^  A  complete, 
fundamental  understanding  of  the  interaction  of  surfaces  with  ion  beam 
radiation,  however,  is  not  yet  available  to  the  materials  scientist  or  process 
engineer. 3 

Surfaces  are  roughened,  not  smoothened,  by  IBE  at  the  sub-micron  level. 

The  influence  of  the  stochastic  nature  of  sputter  deposition  has  been  modeled  by 
Karunsari,  Bruinsma,  aqd  Rudnick"!  by  including  a  shot  noise  term  in  a  linear 
response  theory.  This  model  has  recently  been  modified  for  IBE  by  changing  the 
sign  of  the  derivative  of  the  surface  height  with  respect  to  time,  and  was  found  to 
qualitatively  agree  with  experimental  data.^  In  particular,  the  surface 
roughness  was  found  to  depend  on  the  sample  temperature,  the  ion  fluence.  and 
the  ion  flux  during  IBE. 

The  poresent  paper  is  a  continuation  of  our  prior  work  utilizing  scanning 
tunneling  microscopy  (STM)  to  study  the  surface  morphology  of  ion-etched 
graphite  surfaces. 5-6  In  addition  to  STM.  atomic  force  microscopy  (AFM)  and 
Fourier  transform  infi-ared  (FTIR)  spectroscopy  were  also  used  to  characterize 
the  ion-  etched  surfaces  in  this  work.  The  STM  and  AFM  were  used  because  these 
two  instruments  directly  produce  quantitative  surface  topography  data,  in 
contrast  to  scanning  electron  microscope  and  optical  characterizations.  The 
height  autocovariance  function  for  these  ion-etched  surfaces  was  then  computed  — 
from  the  STM  and  AFM  numerical  data  and  compared  to  model  predictions. 

The  graphite  (0001  '  surfaces  studied  in  these  expenments  were  prepared  by  ~ 
cleaving  highly  oriented  pyrolytic  grapliite  (HOPGi  samples;  these  surfaces  were  — 

atomically  flat  over  tens  of  thousands  of  square  angstroms^  prior  to  ion  beam  Codwg 

d/ot 


exposure.  Graphite  is  also  inert  in  air  and  has  a  very  rigid  lattice,  as 
demonstrated  by  its  high  melting  temperature  of  approximately  3800°C.  This 
means  that  diffusion  after  ion  beam  exposure  is  limited;  the  etching  features  are 
'frozen  in'  to  the  surface  and  can  be  scanned  with  the  STM  or  AFM  long  after  the 
ion  beam  exposure. 

After  cleaving,  the  samples  were  scanned  with  STM  (described  elsewhere^j 
to  insure  that  they  were  atomically  flat  and  then  placed  in  an  ultrahigh  vacuum 
chamber  and  pumped  down  to  about  10*^  torr.  They  were  then  exposed  to  Ar"*"  ions 
accelerated  through  a  potential  of  5  keV  and  rastered  over  a  9  mm^  area  of  the 
graphite  surface  at  an  angle  about  45'  off  the  surface  normal.  A  small  positive 
bias  (45  volts)  was  applied  to  the  sample  to  suppress  secondary  electron  emission. 
In  these  experiments,  the  total  ion  fluence  at  the  center  of  the  sample  was  10^® 
ions/cm^  and  the  ion  flux  conditions  were  6.9  x  10^^,  3.5  x  10^“^,  and  6.9  x  10^"^ 
ions/cm^/sec. 

After  sputtering,  the  samples  were  removed  from  the  UHV  system  through 
a  nitrogen  filled  glovebox  and  stored  there  when  they  were  not  being  scanned  in 
the  STM  or  AFM.  STM  scans  were  collected  in  a  range  of  sizes  from  250A  x  250  A 
up  to  2400A  X  2400A;  AFM  scans  were  collected  at  sizes  from  5000A  x  5000A  up  to 
lOO.OOOA  X  lOO.OOOA.  The  STM  was  operated  in  the  constant  current  or 
topographic  mode  at  a  tunnel  current  of  0.5  nA,  sample  to  tip  bias  of  -100  mV,  and 
in  air  at  room  temperature.  The  AFM  used  was  a  Digital  Instrument.--  Nanoscope 
II  also  operated  in  ambient  air.  Since  oxidation  of  exposed  ion  etched  graphite 
surfaces  has  been  observed  in  samples  exposed  to  ambient  air  for  several  days.^ 
we  compared  the  FTIR  reflectance  spectrum^  (Analect  FTIH  model  fX.'\-530  i  of  an 
IBE  sample  ^re^hly  removed  from  vacuum  to  the  same  sample  after  8  hours  of 
exposure  to  ambient  air.  We  did  not  observe  any  chemical  clianges  at  the  surface, 
and  since  samples  scanned  with  both  STM  and  AFM  were  not  exposed  to  ambient 


3 


air  for  more  than  8  hours,  contributions  to  surface  morphology  due  to  oxidation 
were  neglected. 

In  prior  STM  studies  of  rough  surfaces,  topographs  were  collected  that 
were  shown  to  depend  on  the  shape  and  sharpness  of  the  tips  used  to  collect  the 
data. In  particular,  surfaces  with  sharp  protrusions  can  yield  images  of  the 
sides  of  the  tip,  and  blunt  tips  tend  to  produce  images  which  appear  smoother 
than  they  actually  are  because  of  the  inability  of  the  blunt  tip  to  properly  sample 
small,  sharp  features.  We  have  been  particularly  careful  to  avoid  tip  artifacts  in 
this  study.  In  both  instruments,  topographs  were  collected  with  more  than  one  tip 
and  over  several  areas  of  the  same  sample.  Furthermore,  we  did  not  filter  or 
smoothen  the  data  collected  in  this  study  prior  to  analysis.  Most  of  the  data 
collected  was  reproducible  from  sample  to  sample,  even  with  different  tips. 
Topographs  with  tip  artifacts  were  observed  but  not  analyzed  for  this  study. 

In  these  experiments,  we  varied  the  Ar+  ion  flux  and  compared  the 
topography  of  the  samples  etched  with  10^9  ions/cm^  to  prior  data  at  the  same 
flu.xes  but  at  different  fluences.^  With  the  addition  of  the  AFM,  we  were  able  to 
scan  much  larger  areas  of  the  sputtered  surfaces.  These  extended  range  scans 
turned  out  to  be  crucial  to  propei  lx  identify  and  quantify  the  characteristic  lengths 
and  heights  of  the  fealui'es  pioduced  with  a  fluence  of  10^^  ions'cm-.  Three  large 
scale  AFM  topographs  are  shctwn  in  Fig.  1.  The  average  feature  heights  are 
larger  for  the  higher  two  fluxes,  whereas  the  average  feature  to  feature  length  on 
the  surface  increases  going  from  the  lowest  to  the  intermediate  flux  and  then 
decreases  slightly  at  the  highest  flux.  There  is  also  a  definite  anisotropy  in  the 
average  feature  lengths  in  all  three  topographs  in  Fig.  1.  This  was  verified  to  be 
consistent  over  large  areas  of  the  sample,  and  the  longer  features  were  parallel  to 
the  projection  of  the  incident  ion  beam  on  the  surface.  This  indicates  that 
shadowing'^  played  a  role  in  the  development  of  the  surface  features  at  the  high 
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fluence  used  in  this  study.  No  significant  anisotropy  was  observed  at  the  lower 
fluences  studied  previously,^  which  indicates  that  the  surface  features  must 
reach  some  critical  height  for  a  beam  incidence  angle  of  45°  before  shadowing 
becomes  important. 

We  calculated  the  autocovariance  function  from  the  topographic  data  as 
follows; 


G(lr2-rjl)  =  <h(r2)h(rj)>  -  <h(r)>2  ,  (1) 

where  h(r-)  is  the  height  of  the  surface  at  position  r-  and  L  =  1 1*2  -  I  is  a 

separation  length  on  the  surface.  This  function  provided  a  quantitative 
determination  of  the  base  to  peak  height  variation  and  of  the  short  range  lateral 
correlation  of  the  surface  features.  In  particular,  5=[G(0)]1‘^2  ig  the  interface  width 
or  standard  deviation  of  heights  on  the  surface.  Also,  if  the  initial  drop  off  in  the 
autocovariance  function  is  approximated  by  a  gaussian  distribution,  then  the 
standard  deviation  of  that  gaussian  is  the  characteristic  feature  separation  length 
/.  parallel  to  the  surface  (see  Fig,  2>.  Foi-  the  fluence  of  10^9  ions'cm^,  these  /. 
values  for  the  IBE  surfaces  were  independent  of  whether  they  were  obtained  from 
the  AFM  or  the  STM  topographic  data. 

Table  1  summarizes  these  values  for  the  additional  high  fluence  data 
reported  here  and  for  the  prior  data  at  lower  fluences.^  For  the  entire  set  of  four 
fluences  and  three  fluxes,  the  correlation  lengths  and  interface  widths  clearly 
increase  with  increasing  fluence  for  all  the  fluxes,  which  is  not  surprising.  Since 
the  surfaces  in  all  of  these  experiments  were  initially  atomically  flat,  the  increase 
in  the  interface  widths  is  somehow  related  to  the  amount  of  material  removed 
during  the  sputtering  process.  A  TRIM Monte  Carlo  calculation  predicted  that 
on  average  1.5  C  atoms  are  removed  from  graphite  for  every  incident  Ar'*'  ion 
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under  these  sputtering  conditions.  Therefore,  a  fluence  of  10^^  ions/cm^  would 
correspond  to  the  removal  of  about  4000  monolayers  of  graphite. 

Figui'e  3  shows  the  dependence  of  the  interface  width  as  a  function  of  ion 
fluence  for  all  three  fluxes;  it  is  almost  a  first  power  relationship.  This  result  is 
very  surprising,  because- a  stochastic  model  of  ion  etching  predicts^  that  the 
interface  width,  5,  should  scale  as  the  square  root  of  the  amount  of  material 
removed,  which  is  usually  assumed  to  be  directly  proportional  to  the  ion  fluence; 

5  =  d  (  Q  Y/p  )  0.5  (2) 

where  d  is  the  interlayer  spacing  for  graphite,  Y  is  the  sputter  \ield,  p  is  the  area 
number  density  of  a  graphite  plane,  and  Q  is  the  ion  beam  fluence.  Equation  2  is 
also  plotted  along  with  the  experimental  data  in  Fig.  3  for  Y=1.5  (dashed  line);  it 
clearly  diverges  from  the  data  at  high  fluences.  Other  surface  effects  such  as 
diffusion  and  sputtcn-  redeposition  are  healing  mechanisms;  they  tend  to 
smoothen  features  on  the  surface  and  therefore  could  not  be  responsible  for  an 
increase  in  the  interface  width. 

One  explanation  for  the  data  in  Fig.  3  is  that  the  sputter  yield  is  actually  a 
function  of  the  ion  fluence: 

Y  ~  Yq  (  1  +  (3  Q  )  and  (3) 

6  =  d  [  Yu(  1  +  [3Q)Q 'p]0-5,  (4) 

where  (3  is  a  system-dependent  parameter.  Equation  4  is  also  plotted  with  the 
data  in  Fig.  3  (solid  line'.  The  two  parameters  Y^  and  [3  were  best  fit  to  the 
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experimental  5  vs.  Q  data  (ignoring  the  single  data  point  for  Q  =  IQl®  ions/cm2 
and  J  =  ions/cm^/sec,  which  has  been  previously  shown  to  be  anomalous^j; 

the  values  of  Yq  and  (5  found  in  the  fit  shown  in  Fig.  3  are  1.25  atoms/ion  and  2  x 
lO'lS  cm^/ion  respectively.  This  best  fit  value  of  Yq  is  interesting  because  it  agrees 
fairly  well  with  the  value  of  1.5  independently  calculated  from  the  Monte  Carlo 
simulation,  which  essentially  models  the  sputter  yield  from  a  pristine  surface. 

The  dependence  of  the  sputter  yield  on  fluence  could  physically  correspond 
to  chemical  damage  to  the  graphite  surface.  It  is  important  to  remember  that  the 
surface  before  IBE  in  all  of  these  experiments  was  essentially  a  perfectly  fiat, 
regular  graphite  crystal  layer.  As  the  ions  break  up  the  rigid  pi-orbital  bond 
network  in  the  surface  layers,  the  damaged  surface  could  become  susceptible  to 
enhanced  sputter  removal  by  other  incoming  ions.  This  effect  by  itself,  however, 
would  saturate  at  high  enough  fluences  because  the  depth  of  damaged  layers 
should  tend  to  reach  a  steady  state  value.  However  at  a  fiuence  of  10^^  ions'cm-, 
which  corresponds  to  an  average  of  about  2500  Ar'*'  ions  per  surface  carbon  atom, 
we  see  no  drop  off  in  the  dependence  of  the  interface  width  towards  a  square  root 
dependence.  If  Eq.  4  is  valid,  then  the  sputter  yield  of  a  surface  exposed  to  a 
fiuence  of  10^^  ions 'em-  would  be  about  26  ejected  carbon  atoms  per  incident  Ar'*' 
ion  from  the  fit  shown  in  Fig.  3,  which  is  very  high  compared  to  Yq.  This  suggests 
that  clusters  may  be  ejected  from  heavily  bombarded  surfaces,  and  that  the  total 
amount  of  matei'ial  removed  by  the  10^9  jons/cm^  etch  would  be  closer  to  37.000 
than  to  the  4,000  monolayers  predicted  for  a  constant  sputter  yield.  Since  this  is  a 
substantial  difference  in  total  amount  of  material  removed  by  IBE  at  high 
fluences.  it  is  important  to  measui’e  the  actual  removal  rate  of  material  as  a 
function  of  fiuence.  Such  experiments  are  substantially  more  difficult  to  perform, 
but  we  ai-e  planning  to  examine  this  effect  and  to  also  analyze  other  material 


systems  to  ascertain  whether  or  not  this  is  a  property  unique  to  graphite  or 
general  for  other  materials. 

In  summary,  the  autocovariance  function  has  been  used  to  compress  the 
real  space  height  information  from  STM  and  AFM  topographs  of  ion  sputtered 
surfaces  to  obtain  the  interface  width  and  correlation  length  for  ion  beam  eroded 
surfaces.  We  found  that,  for  a  given  flux,  both  the  correlation  length  and  interface 
width  increase  with  increasing  ion  beam  fluence.  For  a  given  fluence,  however, 
the  interface  widths  appear  to  have  a  non-monotonic  dependence  on  the  ion  beam 
flux,  which  we  do  not  understand  at  present.  Furthermore,  the  interface  widths 
at  a  given  flux  do  not  scale  with  the  square  root  of  the  fluence,  as  is  expected  for  a 
stochastic  process.  This  suggests  that  the  sputter  yield  of  graphite  may  itself  be  a 
function  of  the  fluence.  which  could  result  from  disruption  of  the  graphite  layers 
that  enhances  the  sputter  yield  for  latei'  IBE.  This  apparant  dependence  of  the 
sputter  rield  on  ion  beam  fluence  persists  up  to  fluences  of  10^^  lons  cm^. 
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Figure  Captions 


Figure  1:  Three  2.4ii  x  2.4|a  atomic  force  microscope  topographs  of  a  graphite 
surface  exposed  to  a  total  5  kev  Ar"*"  ion  fluence  of  10 ions/cm^  at  fluxes  of  (ai 
7x1013,  (b)  1.4x10^'^,  and  (c)  7x10^"^  ions/cm^/sec,  respectively.  All  three 
topogi’aphs  are  plotted  so  that  the  x,  y  and  z  axes  are  identically  scaled.  Note  both 
the  asymmetry  of  feature  lengths  and  the  variation  of  average  feature  heights. 

Figure  2:  Plots  of  the  autocovariance  as  a  function  (see  Eq.  1)  of  feature  separation 
length  for  the  three  topographs  shown  in  Fig.  1.  The  interface  width  is  the  square 
root  of  the  y-intercept  (G(0i)  and  the  average  correlation  length  is  the  gaussian 
standard  deviation  of  the  aulocovariance  function  at  low  separation  length  <R). 
The  dashed  curves  correspond  lo  tlie  aulocovariance  parallel  to  the  incoming  ion 
flux  projected  on  the  surface;  the  solid  line  is  the  aulocovariance  perpendicular  to 
the  ion  flux  projection  on  the  surface.  At  the  fluences  examined  earlier,-^  the  two 
autocovariances  did  not  differ  significantly  from  each  other,  but  at  10^^  ions'cm^ 
the  difference  is  pronounced. 

Figure  3:  Log-log  plot  of  the  interlace  width  6  as  a  function  of  the  fluence  Q  for 

the  experimental  data  presented  in  Table  1.  This  plot  also  shows  the  predictions  of 
both  eqns.  2  and  4;  the  deviation  of  the  experimental  data  from  a  0.5  power 
dependence  is  clearly  evident  at  fluences  above  ioix'cm^.  For  the  plot  of  eqn.  4. 
the  best  fit  parameters  Yy  and  [3  were  1.25  iorx'atom  and  2  x  10'^®  cm^/ion. 
respectively. 
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Table  1;  Interface  Widths  and  Correlation  Lengths  (in  A >  as  a  Function  of 
Sputtering  Fluence  and  Flux. 


Fluence 

(ion  s/cm  2) 

Flux  (ions/cm^/sec ) 

Qj  (IOI6) 

Q2  (ioi'7) 

Q3  (lOlSj 

Q4  (10l9j 

Jj  (7  X  lO^^j 

Interface  Width,  5 

1.0 

3.9 

43 

258 

Correlation  Length: 

/.(parallel)^ 

S3 

135 

235 

425 

/.  (perpendiculai'i^^ 

33 

135 

235 

625 

J2  1I.4  X  IOI-*  J 

Interface  Width,  6 

3.3 

8.2 

92 

548 

Correlation  Length: 

/.  (parallel  > 

4t) 

100 

210 

1250 

A  ( perpendicular ' 

40 

100 

210 

1550 

J3  (7  X  IQI-^  I 

Interface  Width.  S 

20.2 

13.1 

31.6 

448 

Correlation  Length; 

/.  (perpendicular) 

20 

33 

220 

525 

/.  ( pai'aHel  1 

2) 

33 

220 

850 

^Correlation  length  along  the 

direction  par 

allel  to  the 

ion  beam  incidence. 

^Cori'elation  length  along  the 

dii'ection  per 

■pendicular 

to  the  ion  beam 

incidence. 
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